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Within a phenomenological quasiparticle model, the quark mass and temperature dependence 
of the QCD equation of state is discussed and compared with lattice QCD results. Different approx- 
imations for the quasiparticle dispersion relations are employed, scaling properties of the equation 
Q-T of state with quark mass and deconfinement temperature are investigated and a continuation to 

O ■ asymptotically large temperatures is presented. 

1 Introduction 

> ■ 

■ Within the past years, physicists aimed at revealing the very nature of strongly interacting matter 
' ■ experimentally by ultra-relativistic heavy-ion collisions performed at CERN-SPS and BNL-RHIC (see 
^ . e. g. [U [2]). At the large local energy densities reached during the collision process, a new state of 

, J ! deconfined matter, dubbed quark-gluon plasma (QGP), is thought to be created. Originally, motivated 

O \ by asymptotic freedom in QCD, the QGP was considered as a weakly interacting gas of quarks (q) 
and gluons (g). However, the success of hydrodynamical concepts (supplemented by assuming fast 
thermalization [31 HI [5J [6] and low viscosity [Zl El [9]) in describing experimental results points to the 
necessity of viewing the QGP rather as a strongly coupled system. 

For hydrodynamic considerations, the QCD equation of state (EoS), which is related to the grand 
& \ potential Q, is of paramount interest. In addition, the EoS is of significant importance for the description 
of the dynamics of the early universe or of compact stellar objects. In general, Q depends on parameters 
being specific for QCD like quark masses m q , flavor number Nf or color number N c as well as on external 
conditions described by temperature T and various chemical potentials {fi q }- From the theoretical side, 
much progress has been achieved in first-principle (lattice) QCD evaluations of the grand potential. 
Previously, the latter were technically restricted to too large m q values of 0(T) translating into a pion 
mass of about 770 MeV [TU] at the deconfinement critical temperature T c . Only recently, smaller quark 
mass parameters were considered [IT] pushing the pion mass to about 215 MeV [121 Nonetheless, 
despite the necessary continuum extrapolation of lattice QCD results, also an extrapolation to physical 
quark masses remains an inevitable step towards obtaining reliable results. 

A variety of phenomenological approaches describes fairly well bulk information found in lattice QCD 
calculations by adjusting suitable parameters. Among these, effective quasiparticle models [TU [T51 dHl 
mug [19], PNJL models [2Q112T1I22], (Polyakov-) quark-meson models [221 HI] and a colored bound state 
model [251 126] have to be mentioned. Here, we employ a quasiparticle model (QPM) for the description 
of QCD thermodynamics. The quark mass dependence is directly implemented in the quasiparticle 
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dispersion relations, which allows for a comparison with lattice QCD results employing different quark 
mass values and enables an extrapolation to the physical limit. The quark mass dependence of QCD 
excitations in one-loop approximation was investigated in [27J. Armed by dispersion relations motivated 
from these considerations, the quark mass dependence of Q is also discussed here. 

Soon, even larger local energy densities may be reached at CERN-LHC, which might provide deeper 
insights into the early formation dynamics of our universe. For pure SU(3) gauge theory, lattice QCD 
thermodynamics was recently studied up to temperatures ~ 10 7 T C [28J. At such asymptotically large 
temperatures, analytical attempts based on perturbative means [291 130] account fairly well for available 
lattice QCD results. Here, we discuss the EoS at asymptotically large temperatures within our QPM, 
both, for pure SU(3) and for Nf = 2 + 1. (Nevertheless, for early universe studies also the influence of 
the heavy quark sector becomes important [3~T].) 



2 Quasiparticle dispersion relations 

In the following, we concentrate on matter anti-matter symmetric systems, represented by a net baryon 
density tib = 0. In this case, the QPM rests on the entropy density s = J2i= g , q s i with 



s t {T)=2t t % r dkk 2 fin 



1 + tie 



where d q = NfN c , d g = N 2 — 1, e q = 1, e g = —1 (for details cf. [151 E2])- This ansatz assumes that 
the thermodynamically relevant excitations at momenta k ~ T are transverse gluons and regular quark 
modes. The pressure p = — Q/V (V is the volume) and other related thermodynamic quantities follow 
from integrating s = dp/dT. The quasiparticle dispersion relations entering Eq. (jTJ can be represented 
by ujg = k 2 + Ug t q(T), where U gq (T) denote temperature dependent self-energies. The self-energies 
can phenomenologically be approximated by 

IL = -G 2 T 2 (N C + ^ 
9 6 V 2 

Ilg = m 2 + 2m 2 + am q mf (3) 

with m 2 = (N 2 — 1) G 2 T 2 / (16 N c ). These expressions (cf. [151 [32] f° r ex = 2) are based on one-loop ap- 
proximations to lowest order in m q in the asymptotic momentum region, where Eq. ([3]) is supplemented 
by a bi-linear term am q mf for including nonzero quark masses according to |33j. G 2 , replacing the 
QCD running coupling g 2 , represents an effective coupling strength parametrized by 

^ 2 -ioop)(C(T)), T>T C , 
G 2 (T) = { ■ / mN / _ (4) 



G( 2 -ioop)(C(T c )) + 6(l-j)+a(l-j) , T<T C 



with ((T) = X(T — T s )/T c which approaches the perturbative region at large T in line with the two-loop 
expression of g 2 . In ((T), X can be related to Aqcd while T s regulates G 2 near T c . For T < T c , G 2 
changes drastically its behavior as dictated by lattice QCD results. 

Recently [27], the quark mass dependence of thermal QCD excitations in the one-loop approximation 
was examined in some detail. The gauge invariant expression for the transverse gluon self-energy at 
asymptotic momenta was found as 



n 9 - i G ^ + l|z(^)) 



(5) 

;mN 12 (m a \ 2 r 00 , a 2 1 

l(m T) = — (— i] / da - ; = , . 6 

1 q ' ' 7T 2 V T J Jo ^YTa 1 (1 + e v / T+^m 9 /T) 
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The power expansion of I(m q /T) for small m q /T involves a term resembling chiral logarithms, cf. [27J. 
Neglecting any quark mass dependence by setting m q /T — > 0, one gets T{m q /T) — > 1 and thus Eq. (j2J) 
is reproduced. This approximation is supported to some extent by the one-loop results [27], as the 
energy of transverse gluon excitations increases only by a tiny amount when decreasing m q . 

In the case of regular quark excitations, U q = m 2 + m 2 T q {k,m q ) is found in [27] . where T q (de- 
pending non-trivially on k and m q ) encodes the effects induced by the thermal medium. Straightfor- 
ward evaluation of this expression in the asymptotic momentum region yields Il g = m 2 + 2M^_ with 
= ^m 2 {T{m q /T) + 2) and X{m q /T) given in Eq. Thus, in the limit m q — > 0, reduces 
to m 2 . Obviously and in contrast to Eq. ([3]), a bi-linear term relating m q and rrif is not present in 
this result. However, for momenta k ~ T and small m q , a bi-linear term can be motivated from the 
one-loop approximations. In general, one can represent J-' q {k,m q ) = a{k,m q }T q {k = 0,m q ), where, for 
small m q (and fixed k), a is rather independent of m q and a number between 1 and 2. Approximating 
T q {k = 0,m q ) for small but nonzero m q and small g in line with a generalization of [55] . 

U q = m 2 q + a (--mj + m q M + + M^j (7) 

is obtained for an approximation of the asymptotic quark self-energy. With this ansatz, the thermal 
quark mass effectively decreases with decreasing quark mass parameter, as in the case of Eq. ([3]). In 
fact, both approximations Eq. ([3]) and Eq. ([TJ yield the same expression for the quark self-energy in 
the limit m C 1 and for a = 2 reading IT^ = 2m q rrif + 2m 2 . 

3 Quark mass extrapolation 

Benchmark of the considerations is the scaled entropy density s/T 3 as a function of T/T c for Nf = 2 + 1. 
In [TO], a continuum estimate for s/T 3 is given by reporting lattice QCD results for the energy density e 
and the interaction measure A = (e— 3p), where s/T 3 = (4e — A) / (3T 4 ), for fairly large and temperature 
dependent quark mass parameters m u ^/T = 0.4 and m s /T = 1. Recently [HI [12], s/T 3 was calculated 
for much smaller quark masses. Corresponding to [31] , the quark mass parameters used in [TTj fl~2] can 
be approximated by m u /T = m d /T = A /T 2 + A 1 /T + A 2 with A = 0.298- 10~ 3 GeV 2 , A 1 = 1.664- 10~ 3 
GeV and A 2 = 0.002 for T given in GeV and m s = 10 m u . As s/T 3 in [12] is given as a function of T 
in MeV, we scale by T c = 190 MeV in line with [12], where T c = (190 ± 5) MeV is reported giving rise 
to an estimated error in T/T c according to AT C = 5 MeV. Despite, T c is assumed to be approximately 
quark (and related pion) mass independent (cf. [35] and references therein). 

In the following, the QPM armed by the two different approximations for the dispersion relations 
discussed in Sec. [2] is adjusted to the lattice QCD results from [TP] . The application of Eqs. ([2]) and ([3]) 
in the quasiparticle dispersion relations is denoted by "Fit 1", while the use of Eqs. ©-([Tj) is associated 
with "Fit 2". Then, by extrapolating to the mass set-up employed in [12], the proper implementation 
of the quark mass dependence in the QPM can directly be tested. For simplicity, any conceivable m q 
dependence in G 2 is naively neglected. 

Starting with Fit 1, i. e. Eqs. (T5]) and Q as approximation of the self-energies, we choose a = 1, 
here, as reasonable value (cf. Sec. [2]). The corresponding QPM parameters of G 2 read T s = 0.7 T c , A = 5, 
a = —426, b = 403.3 describing s/T 3 from [10j impressively well as exhibited by the lower dashed curve 
in Fig. [1] (left panel). Extrapolating to smaller quark masses as used in [12], the corresponding lattice 
QCD results are fairly well reproduced, cf. upper dashed curve in Fig. [1] (left panel). Decreasing m q , the 
thermal gluon mass remains unaffected, while the thermal quark mass is reduced according to Eq. 
Thus, the entropy density increases with decreasing quark mass. 

Note that an inclusion of the term am q mf of significant strength seems to be mandatory. Neglecting 
this contribution by setting a = 0, any adjustment of QPM parameters to [10| fails in describing the 



3 




Figure 1: Left: scaled entropy density s/T 3 as a function of T/T c for Nf = 2 + 1. Circles [10] 
and squares [12] exhibit lattice QCD results for different quark mass set-ups. Dashed and 
solid curves depict QPM results for Fit 1 and Fit 2, respectively, employing different ap- 
proximations for the quasiparticle dispersion relations as defined in Sec. [2j Lower curves 
represent adjustments to [TU] while upper curves represent quark mass extrapolations in line 
with [12]. Right: scaled interaction measure A/T 4 = (e — 3p)/T 4 as a function of T/T c . 
Squares denote lattice QCD results from [121 EE], dashed and solid curves corresponding 
QPM results for Fit 1 and Fit 2, respectively. Dotted curve depicts the fuzzy bag model [36] 
result for T > 1.08 T c . 



quark mass extrapolation to |12j . This highlights the crucial role of the bi- linear term for the quark 
mass dependence of the phenomenological quasiparticle dispersion relations. For larger a, G 2 has to 
take a smaller value at fixed T for describing [10J. However, a cannot be chosen arbitrarily large, as a 
dominance of am q mf might lead to an overestimation of the lattice QCD results [12J when extrapolating 
to the according quark mass values. We mention that using, instead, constant quark mass parameters 
m u = md = MeV and m s = 55 — 90 MeV (which might be considered as physical limit) does not 
noticeably change the upper dashed curve on the scale exhibited in the left panel of Fig. [fl 

In the case of Fit 2, i. e. Eqs. ([5]) - ([7]) as ansatz for the self-energies, we first note that the thermal 
gluon mass increases with decreasing m q according to Eq. ([5]), while the behavior of the thermal quark 
mass is, in general, influenced by two counter-acting effects: M + increases with decreasing m q whereas 
the terms solely oc m q decrease. Thus, depending on an appropriate choice for a, U q might decrease 
with decreasing m q mostly as a result of am q M + . Consequently, a larger value for a has to be expected 
compared to Fit 1 in order to saturate the quark mass behavior of s/T 3 from [12J. Here, for example, we 
choose a = 2. The corresponding QPM parameters, being similar to Fit 1, read T s = 0.68 T c , A = 4.77, 
a = —426 and b = 403.3 when adjusting to s/T 3 from [TO], cf. lower solid curve in Fig. [1] (left panel). 
The extrapolation to smaller quark masses in line with [12] is exhibited by the upper solid curve in the 
left panel of Fig. [lj In both cases, Fit 1 and Fit 2 (on the given scale solid and dashed curves lie almost 
on top of each other in Fig. [I]- left panel), agreement with the corresponding lattice QCD results is 
found. 

In addition, in [121 IS]; the scaled interaction measure A/T 4 was calculated. The QPM results, 
according to Fit 1 and Fit 2 for the small quark mass set-up, are depicted by dashed and solid curves 
in the right panel of Fig. (H respectively, finding an overall good agreement. The necessary pressure 
integration constant [15], denoted by B(T C ), reads B(T C ) = 0.65 T 4 . Nonetheless, the QPM m q extrap- 
olations overshoot somewhat the lattice QCD results around T c (Fit 1 more than Fit 2) and exhibit 
small deviations also at larger temperatures. 
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Contrary, in a fuzzy bag model approach [36], lattice QCD results of A/T 4 for T > 1.2 T c can 
perfectly be described, cf . dotted curve in Fig. [TJ (right panel) . The fuzzy bag picture originally accounts 
accurately for the plateau observed in the scaled interaction measure A/T 2 as a function of T (cf. dotted 
curve in Fig. [2]- left panel), which is parametrized by A/T 2 = 2Bf + ABmit/T 2 . Considering leading 
non-perturbative contributions to the pressure to be given by a temperature dependent (fuzzy) bag 
constant Bf mimicking a gradual rather than an abrupt transition from the confined phase to a nearly 
perturbative phase, this picture represents a generalization of the MIT bag model. The fit parameters 
reproducing A/T 4 for T > 1.2 T c from [El [13] in Fig. □ read B f = 0.135 GeV 2 , B MIT = 0.0009 GeV 4 . 
Nonetheless, this approach neither describes the behavior of the interaction measure in the transition 
region nor incorporates explicitly quark mass effects. 



4 Scaling properties of the equation of state 

Discussing the scaling properties of the QPM EoS with m q and the value of T c , one first notes that the 
pressure reformulated in a dimensionless fashion reads 



m = = E e vH dxx 2 ln(l + e ie 



i=g,q 



JO 



2tt 2 £Wi V 1 dr)J Wo ./^5 e V^ + 





with IT = n,(T)/T 2 , = n^TO/T' 2 , real number b = B{T C )/T* and £ = T/T c , rj = T'/T c . Eq. © 
depends explicitly on m q via IT. For m q = eT as used in JTU], IT (TT^) is a function of £ (r/) only 
because the entering G 2 depends on T c only via £ (rj) (note that AT S /T C is a parameter in the QPM). 
Thus, $(£) displayed as a function of £ is independent of the explicit value of T c . Consequently, 
related thermodynamic quantities like s/T 3 and A/T 4 show the same independence of T c . Nonetheless, 
approximating m q as advocated above in line with [HI [12], Il q explicitly depends on T, and thus, $(£) 
does not exhibit the discussed T c independence. However, numerically even a variation in T c by 100 
MeV turns out to imply negligible effects on p/T 4 , s/T 3 or A/T 4 . 

The interaction measure scaled by T 2 and depicted as a function of T shows a T c dependence, 
cf. Fig. [2] (left panel). A shift in T c by 5 (—5) MeV in line with [12] results in a shift of the curve A/T 2 
to the right (left) in the transition region. In addition, the curves suffer a parallel shift up (down) by 
about 6-7% for larger temperatures. Note that the deviations in A/T 4 observed in Fig. [TJ between QPM 
and [T2"| [T3"] for larger T are transparently quantified by deviations from the plateau behavior in Fig. [2] 
(left panel). 

The parallel shift in A/T 2 can qualitatively be understood from the scaling behavior of the EoS with 
T c . Reformulating p = T 4 £($£ 3 ) and imposing $(£) = $(£')> cf. Eq. ([8]), p and thus e = — p + £<9p/<9£ 
change when changing T c to T' c (i. e. £ to £') according to p' = p{T' c /T c ) A and e' = e(T^/T c ) 4 . For 
the EoS in the form p(e), this implies that the linear section at larger e, which can be approximated 
by p(e) = ae + po, becomes p'(e') = ae' + (T^/T c ) 4 p with the same slope a but different off-set 
p' = (T c '/T c ) 4 p - For T' c > T c , the linear section of p(e) is, thus, parallely shifted downward whereas for 
T' c < T c it is shifted upward. A similar behavior is observed for the EoS in the form p/e as exhibited in 
Fig. [2] (right panel). 

Other regions of p(e) might need to be approximated differently, say for instance, by p(e) = ay/e+po 
close to the transition region. Changing again T c to T c ', the EoS changes into p'(e') = a'^fd + p~o with 
a' = (T c /Tc) 4 a and po = {T' c /T c ) A p~o implying a change in the off-set p~o but also a flatter curve, a' < a, 
for T' c > T c and a steeper curve, a' > a, for T' c < T c . This leads, now for p/e, to a change in the order 
of curves as evident from Fig. [2] (right panel) at lower energy densities. 
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Figure 2: Left: scaled interaction measure (A/T 2 )/GeV 2 as a function of T/GeV. Squares 
exhibit lattice QCD results from [121 [13]. Dashed, solid and dash-dotted curves depict QPM 
results according to Fit 2 using T c = 185, 190, 195 MeV, respectively. Dotted curve shows 
the directly fitted fuzzy bag model [35] result. Right: EoS in the form p/e as a function of 
e/(GeV/ fm 3 ) and its T c dependence. Lattice QCD results (squares) from [37], curves depict 
QPM results with the same line code as in the left panel. 



Likewise, one may discuss the m q dependence of the EoS. As empirically evident from Fig. [Q the 
m q dependence of s/T 3 can be parametrized by s'/T 3 = m(C,)s/T 3 with T dependent function m. 
Similarly, p and e (as discussed in |32j) behave according to p' = m{^)p and e' = m(£)e + £p<9m(£)/<9£ 
with m q . Numerically, one finds m(£) = 1.107... 1.081 in the interval £ = 1.5. ..3, i. e. small changes 
in m(£) with £, whereas for small £ a larger m(0.8) = 1.33 and m(0.9) = 1.51 is found. Thus, 
m(£) varies more sizeably with £ for small £. Approximating p(e) again by a linear function, p' = 
ae' / (l + 7^ d^(£) j _|_p Qm/ / _|_ g£ ) is found. Consequently, changes in p(e) with m q are negligible 
for larger £, where m(£) ~ 1 is almost constant, whereas they are mostly visible in a region of £, where 
m(£) varies most rapidly, i. e. for £ < 1. 



5 Equation of state at asymptotic temperatures 

Recently [28], the pressure for pure SU(3) gauge theory became available in a large temperature interval 
between T c and 10 7 T C confirming old results [38J between T c and 5T C . As p changes most rapidly in 
the transition region, we fit the QPM parameters to [38] and then continue to larger T. The result 
for p/p SB , where p denotes the Stefan-Boltzmann pressure, is exhibited by the solid curve in Fig. [3] 
(left panel). While the QPM describes appropriately the regions around T c and at asymptotically large 
temperatures, say for T > 200 T c , it underestimates [28] in the intermediate T region between 3 T c 
and 200 T c . Although one might find an exotic QPM parametrization for intermediate and large T 
reproducing [28] and delivering similar results as perturbative QCD [303 SI] > cf • dashed curve in Fig. [3] 
(left panel), such a description fails for T < 2 T c . This discrepancy might be viewed as a hint for missing 
degrees of freedom in the intermediate T region. Considering a hard-thermal-loop inspired extension of 
the QPM [ID], including in addition longitudinal gluon (plasmon) and Landau damping contributions, 
the situation observed in the left panel of Fig. [3] is not measurably improved. This is due to the fact 
that the effect of a negative partial plasmon pressure decreasing for T — > T c + is compensated by the 
interwoven increasing partial transverse gluon pressure (driven by the Landau damping) for T — > T c + , 
while the plasmon contribution vanishes for increasing T. Nonetheless, to resolve this issue, additional 
lattice QCD studies would be desirable. 
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Figure 3: Left: comparison of QPM (solid and dashed curves) with lattice QCD results 
(circles [38] and squares [28]) for p/p SB as a function of T/T c > 1 of a pure gluon gas. The 
QPM is adjusted to [38] with T s = 0.855 T c , A = 7, B{T C ) = -0.03 T c 4 (solid curve) or to [28] 
with T s = 0.97T C , A = 150, B(T C ) = 5T C 4 (dashed curve), where T c = 271 MeV [39J . Dashed 
curve is similar to perturbative QCD results reported in [3UI [31]. Right: scaled interaction 
measure (Fit 2 in Fig. [TJ right panel) as a function of T/T c continued to larger T. 



Turning to Nf = 2 + 1 at asymptotic temperatures, the continuation of the scaled interaction 
measure (cf. Fit 2 in Fig. [T]- right panel) to larger T is exhibited in Fig. [3] (right panel). For large T, 
A/T 4 approaches zero logarithmically in line with the temperature dependence of G 2 . While already 
for T ~ 10 T c the conformal limit e = 3p is approximately reached, pjT A still exhibits deviations from 
p /T A of about 10% at T = 50 T c and of 4.5% at T = 10 5 T C . As the dynamically generated thermal 
gluon and quark masses exhibit a behavior ~ GT in the QPM, they are of the order of 12 TeV and 8 
TeV, respectively, at T = 10 5 T c , invalidating the naive picture of weakly coupled quarks and gluons 
with negligible masses. 

6 Conclusion 

In summary, we study the quark mass and temperature dependence of the QCD equation of state. 
We utilize a quasiparticle model by employing two different expressions for the quasiparticle dispersion 
relations, which explicitly depend on m q and are based on one- loop QCD approximations [27] . In both 
cases remarkable agreement with first-principle lattice QCD results is achieved when extrapolating in 
the m q parameter space. Scaling properties of the QCD EoS with m q and T c are discussed and the 
EoS is continued to asymptotically large T, where very hot QCD matter may be viewed as composed 
of rather heavy quasiparticle excitations. 

The authors thank E. Laermann and Z. Fodor for valuable discussions. The work is supported by 
BMBF 06DR136 and EU I3HP. 
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